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1 .  INTRODUCTION 


Typically,  the  report  of  an  unusual  seismic  event  in  a  remote  area  will  generate  a  thorough 
investigation  of  previous  historical  seismic  activity  in  that  area.  For  historically  reported  activity, 
error  analysis  is  complicated  because  phases  for  events  occurring  in  different  time  periods  are 
reported  by  different  stations,  arrivals  are  read  by  different  individuals,  and  locations  have  been 
determined  by  different  methods.  For  this  reason  it  is  useful  to  develop  a  framework  for 
estimating  location  uncertainty  that  is  not  dependent  on  a  formal  statistical  analysis  such  as 
proposed  by  Flinn  (1965)  or  Boyd  and  Snoke  (1984). 

The  most  common  methods  for  locating  seismic  events  all  use  least  squares  to  compare 
arrival  time  observations  with  times  calculated  using  a  regional  or  global  earth  model.  In  practice 
when  different  individuals  or  organizations  independently  locate  seismic  events  they  obtain  slightly 
different  results  for  three  essential  reasons: 

1)  Readings  of  the  arrival  times  differ,  due  to  unavoidable  subjective  differences  in 
identifying,  picking,  and  reading  phase  arrivals; 

2)  There  exist  systematic,  station-dependent  effects  on  seismic  travel  times  of  a  magnitude 
which  depends  specifically  on  the  details  of  the  location  method,  which  includes  explicit  or 
implicit  decisions  about  the  relative  weights  to  be  assigned  to  different  arrivals  or  to 
different  observing  stations; 

3)  The  appropriate  velocity  model  for  calculating  travel  times  is  always  imperfectly  known. 

The  present  study  will  evaluate  the  relative  contribution  of  each  of  these  factors  to  locatitm 
uncertainty  by  careful  analysis  of  seismic  locations  in  three  quite  different  geographic  regions 
(Table  1): 

1)  Macquarie  Ridge  -  This  remote  region  exteixls  to  about  1S(X)  km  south  of  New  Zealand  and 
possesses  a  high  rate  of  shallow  earthquake  activity,  but  the  only  available  phase  data  is 
from  teleseismic  stations; 

2)  Bucaramanga  Nest,  Colombia  -  This  is  an  unusual  intermediate-depth  nest  of  earthquakes, 
observable  both  at  regitxial  and  at  teleseismic  stations. 

3)  Efate  Nest,  Vanuatu  -  There  exists  a  peculiar  concentration  of  shallow  earthquake  activity 
just  to  the  west  of  Efate  Island,  Vanuatu  (formerly,  the  New  Hebrides).  To  evaluate  this 
activity  we  have  observations  from  a  regional  seismic  network  and  from  temporarily 
deployed  ocean  botttxn  seismograph  (OBS)  stations. 
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Table  1.  Comparative  Characteristics  of  Three  Hypocentral  Groups 


Geographic  Region 

Focal  Deptii 

Data  Type 

Data  Source 

Macquarie  Ridge 

Shallow  (~20  km) 

Teleseismic 

WWSSN  Microfiche 

Vamiatu 

Shallow  (~30  km) 

Local 

Field  Observations 

Bucaramanga 

InteimediatB  (~160  km) 

Telesdsmic 

Regional 

WWSSN  Microfiche 

Although  both  the  character  of  the  data  and  the  nature  of  the  location  problem  differs  for 
each  of  these  three  data  sets,  in  each  case  there  are  two  basic  types  of  data  available  fen*  analysis; 

1)  Arrival  time  data  -  This  includes  arrival  times  read  by  the  author  and/or  a  close  colleague 
which  we  can  compare  statistically  with  arrival  times  reported  elsewhere; 

2)  Seismic  event  locations  -  These  include  locations  determined  by  the  author  from  specially 
selected  groups  of  stations  for  cmnpariscxi  witii  other  reported  locations. 

For  all  newly  determined  locations  we  use  the  TexFlex  program  developed  specifically  for 
titis  study  and  described  in  detail  elsewhere  (Brohlich,  1992;  1993).  One  objective  of  this  study  is 
to  test  the  capability  of  the  TexFlex  program  using  a  broad  variety  of  input  data.  The  TexFlex 
program  is  capable  of  locating  seismic  events  using  either  user-supplied  flat-eartii  velocity  models, 
using  teleseismic  look-up  tables  like  the  Jeffreys-Bullen  (1970)  tables,  or  using  the  IASPEI91 
teleseismic  software  described  by  Kermett  (1991)  and  Kennett  and  Engdahl  (1991).  TexFlex  can 
perform  either  single-event,  joint  hypocenter  determination  (JHD),  or  joint  epicenter  determination 
(JED)  relocations.  The  program  allows  the  user  to  choose  a  variety  of  phase  weighting  and 
station-grouping  schemes.  Individuals  interested  in  obtaining  the  current  version  of  the  TexFlex 
software  and  an  explanatory  manual  should  contact  the  author  by  email. 

The  next  section  will  summarize  the  arrival  time  data  and  locations  for  each  of  the  three 
geographic  regions.  Then,  the  sections  following  will  make  some  comparisons  between  location 
uncertainties  for  seismic  events  in  the  various  regions,  and  draw  some  general  conclusions  about 
strategy  for  locating  seismic  events  using  travel-time-based  methods.  Finally,  the  last  section  will 
summarize  some  implications  of  the  location  analysis  in  each  geographic  region  which  may  be 
impmtant  for  individuals  interested  in  the  specifle  ^graphic  areas  studied. 


2 


2.  ANALYSIS  OF  INDIVIDUAL  EARTHQUAKE  GROUPS 
A.  Macquarie  Ridge 

Background:  The  geographic  fegicm  between  New  Zealand  and  about  60"S  is  lenowiKd 
for  its  remoteness,  its  poor  weather,  and  to  seismologists  (Ruff  et  al.,  1989)  for  its  abundance  of 
large  shallow  earthquakes  (Figure  1).  Most  of  the  seismic  activity  falls  along  a  shallow 
bathymetric  feature  called  the  Macquarie  Ridge  which  is  the  boundary  between  the  Pacific  and 
Australian  plates. 

The  Macquarie  Ridge  is  unusual  seismically  in  that  it  possesses  earthquakes  having  a  broad 
variety  of  focal  mechanisms,  including  strike  slip,  normal,  and  thrust  mechanisms,  as  well  as 
events  with  distinctly  non-double-couple  mechanisms  (Frohlich  et  al.,  1989).  However,  the  vast 
majority  of  all  earthquakes  south  of  about  48**$  have  mechanisms  with  approxinuttely  hcnizontal  P 
axes  and  approximately  vertical  B  axes,  consistent  with  predominantly  strike-slip  motion  along  die 
Macquarie  Ridge  (Figure  2) 

There  has  been  special  interest  in  this  region  (e.g.,  Ruff,  1990)  because  on  23  May  1989 
the  largest  (Mw  =  8.2)  submarine  strike-slip  earthquake  in  the  wodd  this  century  occurred  there  at 
about  52®S,  160®E.  Recently  Dm  (1993)  used  phase  arrivals  reported  to  the  International 
Seismological  Centre  (ISC)  and  relocated  earthquakes  along  the  Macquarie  Ridge  using  a  JED 
method.  Previously,  the  locations  reported  by  the  ISC  (Figure  3)  were  the  best  available  for 
seismic  events  in  this  area  occurring  since  1964.  hi  the  next  section  we  shall  explicitly  compare  the 
locations  reported  by  the  ISC  and  those  reported  by  Das  (1993)  to  locations  determined  explicitly 
for  this  study  using  a  select  group  of  statiems. 

Data  Analysis:  Because  the  Macquarie  region  is  so  remote  from  most  of  the  world's 
seismic  stations,  only  the  largest  earthquakes  are  recorded  by  most  stations  in  the  global  network, 
and  most  earthquakes  are  located  using  only  a  few  regitxial  stations.  For  this  reason  we  determiiie 
relocations  using  only  a  small  group  of  stations,  selected  individually  because  they  are  as  sensitive 
and  reliable  as  any  available  stations,  and  selected  as  a  group  to  represent  the  best  possible 
azimuthal  covoa^  for  events  in  the  Macquarie  region. 

The  final  station  set  chosen  included  17  stations  (Figure  4).  This  included  two  Canadian 
stations  (FCC  and  MBC)  at  the  appropriate  distance  to  record  PKP  phases.  However,  because 
these  station  consistently  had  large  residuals,  and  because  phases  observed  there  were  often  quite 
weak  for  Macquarie  events,  FCC  and  MBC  were  ntx  used  for  the  final  relocations. 
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Figure  1.  Map  of  locations  of  earthquakes  with  magnitudes  excee^ng  6.8  reported  since  1904 
a^g  die  Macquarie  Ridge.  Fot  earthquakes  occurring  since  1S>77  locations  and  focal  mechanisms 
are  mnn  the  Harvard  CMT  catalog,  and  magnitudes  are  moment  magnitudes  Mw  Ptior  to  1977 
locations  and  mechanisms  are  as  reported  by  Ruff  et  al.  (1989)  if  available.  All  remaining 
magnitudes  are  mg  finom  Abe  (1981)  remaining  locations  are  from  Abe  (1981). 


4 


i  I  I  I  i _ I _ I — I — I — I — I — I — I 

156  160  164  168 


Figure  2.  Map  of  locations  of  earthquakes  with  CMT  reponed  by  Harvard,  and  composite 
mechanisms  for  five  geographic  groups  altmg  the  Macquarie  Ridge.  On  focal  spheres,  triangles 
are  T  axes,  octagons  are  P  axes,  and  pluses  are  B  axes  for  CMT  of  earthquakes  in  adjacent 
geographic  region  surrounded  by  box,  while  the  plotted  focal  mechanism  remsents  a  composite 
mechanism  determined  by  adding  for  the  individual  earthquakes.  For  CMT  with  significant 
non-double-couple  compcments  (ratio  of  minimiun  to  maximum  principal  moments  exceeding 
0.15),  only  the  well-determined  polar  axis  of  the  CMT  is  plott^.  Except  in  the  Fiordland 
subregion,  note  that  die  predominant  mechanism  type  is  strike-slip,  with  approximately  horizontal 
P  axis  and  qiproximately  vertical  B  axis. 
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Figure  3.  Map  of  earthquake  locations  reported  by  the  ISC  along  the  Macquarie  Ridge 
earthquakes  occulting  after  1964  and  having  body-wave  magnitude  mb  exceeding  5.0. 
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Using  the  WWSSN  microfilm  collection  at  Columbia  University,  we  obtained  copies  of  P 
arrivals  at  tlwse  stations  for  as  many  earthquakes  as  possible.  However,  in  practice  it  was  only 
possible  to  find  readable  arrivals  at  the  majority  of  stations  for  earthquakes  with  magnitudes 
exceeding  mb  of  5.3  or  so.  Altogether  we  obtained  phases  for  S3  earthquakes,  including  36  with  a 
sufficwnt  data  quality  diat  we  felt  we  could  relocate  them  more  accurately  than  did  the  ISC  (Figure 
Sand  Table  2). 

Of  these  36  earthquakes,  '  were  also  relocated  by  Das  (1993),  who  used  the  same  phase 
arrival  data  as  reported  by  the  ISC^  but  obtained  different  locations  because  she  determined  station 
corrections  using  the  JED  method  of  Dewey  (1971).  A  comparison  of  the  locaticms  reported  by  the 
ISC,  by  Das  (1993)  and  this  study  (Figure  6)  shows  that  median  discrepancy  between  the  locations 
of  the  ISC  and  this  study  was  Disc-JHD  =  23  km,  while  the  median  discrepancy  between  the  ISC 
and  Das  (1993)  locations  was  Disc-Das  =  12  km.  In  a  few  cases,  the  di^erence  in  location  was  80 
km  or  more. 

For  the  arrivals  for  S3  earthquakes  reread  especially  for  this  study,  37%  agreed  with  the 
ISC  reported  arrivals  within  0.S  seconds  (Table  3).  However,  17%  differed  from  the  ISC  reported 
arrival  by  1.0  seconds  or  more.  The  mean  difference  between  phase  arrivals  read  by  us  and 
reported  by  the  ISC  was  0.37  seconds.  Moreover,  15%  of  the  phases  we  found  and  reread  were 
not  reported  to  the  ISC.  Also,  there  were  ISC  reported  phase  arrivals  for  24%  of  the  phases  which 
we  were  unable  to  find  after  a  careful  search  of  the  WWSSN  microfiche  seismograms. 

Using  the  phases  reread  firom  the  WWSSN  seismograms,  we  relocated  events  in  both  the 
single-event  and  JED  mode,  and  using  both  IASPEI91  and  Jeffreys-Bullen  (1970)  travel  time 
tables.  For  IASPEI91  travel  times  the  use  of  JED  reduced  the  variance  of  travel-time  residuals 
(Table  4)  from  1.25  sec^  to  0.83  sec^,  corresponding  to  a  reduction  in  variance  of  34%.  For 
Jeffreys-Bullen  travel  times,  the  single-event  variance  was  slightly  larger  (1.35  sec^),  but  the 
variance  for  JED  was  somewhat  smaller  (0.71  sec^),  for  a  reduction  of  48%. 

B .  Bucaramanga,  Colombia 

Background:  About  160  km  beneath  Bucaramanga,  Colombia  (Figure  7),  there  exists  a 
remarkable  concentration  of  earthquake  activity  known  as  the  Bucaramanga  nest  (Trygvasson  and 
Lawson,  1970;  Dewey,  1972).  Since  1964  the  ISC  has  reported  an  average  of  about  20 
earthquakes/year  with  magnitude  mb  of  4.7  or  above  occurring  within  the  nest.  Yet,  careful 
studies  using  data  from  a  temporary  local  network  suggest  that  the  active  area  has  dimensions  of 
only  about  4-8  km  and  a  total  volume  of  about  1 1  km^  (Pennington  et  al.,  1979;  Schneider  et  al., 
1987;  Frohlich  and  Kadinsky-Cade,  1993). 
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Table  2.  Epicenters  near  the  Macquarie  Ridge  determined  by  using  the  JED 
method  and  arrival  times  read  personally  by  the  author. 


Dale 

Hr  Min  Sec 

1  jititiiH#! 

®s 

Longitude 

®E 

Dqah 

km 

rms 

sec 

11/8/66 

2  43  56.01 

49.109 

163.905 

20.00 

1.20 

4/5/66 

1157  41.83 

54.872 

158.610 

20.00 

1.32 

4/28/66 

1  15  38.44 

49.148 

164.161 

20.00 

1.04 

S/25/66 

13  2057.00 

52.677 

160.362 

20.00 

0.37 

6/28/67 

14  34  5.54 

47.004 

165.408 

20.00 

0.87 

9/20/67 

9  39  15.36 

49.676 

163.653 

20.00 

0.74 

9/20/67 

10  30  55.59 

49.641 

163.514 

20.00 

0.60 

9/23/67 

7  2  7.00 

49.555 

163.834 

20.00 

0.40 

5/8/68 

11  0  7.76 

57.882 

158.004 

20.00 

0.56 

9/25/68 

7  2  51.56 

46.363 

166.553 

20.00 

0.56 

6/11/70 

16  46  37.87 

59.367 

156.511 

20.00 

1.26 

3/23/72 

23  10  0.74 

45.266 

166.400 

20.00 

0.65 

4/1/72 

23  51  24.31 

49.347 

163.778 

20.00 

0.50 

4/2/72 

0  17  45.37 

49.378 

163.876 

20.00 

2.10 

12/24/72 

20  3058.81 

52.259 

160.765 

20.00 

1.52 

6/7/73 

2  43  31.83 

53.890 

158.956 

20.00 

0.43 

10/19/73 

0  13  1.68 

54.701 

158,371 

20.00 

0.66 

4/6/74 

12  16  24.72 

48.814 

164  470 

20.00 

0.83 

6/5/74 

1144  22.12 

54.719 

157.258 

20.00 

0.23 

9/10/75 

15  48  54.96 

49.413 

163.824 

20.00 

0.32 

9/16/75 

0  0  50.26 

47.094 

165.242 

20.00 

0.83 

1/7/76 

6  59  12.32 

47.119 

164.899 

20.00 

0.89 

7/21/77 

11  53  22.21 

53.821 

158.921 

20.00 

0.82 

7/21/77 

12  34  32.95 

53.919 

158.117 

20.00 

0.42 

5/24/78 

6  12  2.48 

52.972 

159.560 

20.00 

0.83 

8/28/78 

8  4  14.06 

48.842 

165.209 

20.00 

0.80 

5/25/81 

5  25  14.61 

48.580 

164.767 

20.00 

1.29 

5/30/81 

9  47  19.52 

49.009 

164.488 

20.00 

0.40 

5/23/84 

5  16  37.86 

51.881 

160.883 

20.00 

0.80 

8/2/65 

13  19  55.73 

56.099 

156.888 

20.00 

1.28 

9/12/64 

22  7  3.74 

48.923 

164.341 

20.00 

0.42 

9/03/87 

640  14.05 

58.803 

158.607 

20.00 

0.49 

9/03/87 

8  1  36.46 

59.452 

159.027 

20.00 

0.55 

2/07/80 

10  59  13.70 

54.194 

158.943 

20.00 

0.46 

4/02/72 

0  39  2.10 

49.407 

163.905 

20.00 

0.53 

6/26/71 

22  2  24.30 

51.510 

167.593 

20.00 

1.63 
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Diff.  in  Longitude  (km) 


Figure  6.  Differences  in  locaticms  for  35  epicentn^  determined  in  this  study  (JHD),  and 
those  reported  by  the  ISC  and  by  Das  (1993). 
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Table  3.  Time  differences  between  phase  arrival  picks  reported  by 
differoit  sources.  Results  are  presented  as  percentages  of  585  iHcks  at  16 
stations  for  53  earthquakes  along  the  Macquarie  Ridge,  of  393  pidcs  at  18 
stations  for  27  earthquakes  near  Bucaramanga,  Colombia,  and  of  5,991 
pidu  at  20  stations  for  earthquakes  near  Efate,  Vanuatu.  For  Bucaramann 
and  Macquarie  Ridge,  sources  of  picks  are  arrivate  reported  by  the  ISC, 
and  arrivals  picked  by  the  author  and  his  colleagues.  For  Efate,  the 
arrivals  were  picked  independently  by  two  different  colleagues. 


Time  Difference 
(sec) 

Macquarie  Ridge 
Frac^ofTo^ 

% 

Bucaramanga 
Fraction  of  Total 
% 

^ate 

Fraction  of  Total 
% 

0.0  -  0.5 

37 

38 

87 

0.5  -  1.0 

17 

13 

10 

1.0  -  2.0 

10 

4 

2 

2.0  -  4.0 

4 

1 

0 

>4.0 

3 

2 

1 

No  ISC  time  repotted 

15 

26 

- 

No  phase  found 
by  author 

24 

16 

• 

Mean  difference 

0.37'^  sec 

0.39^  sec 

0.21*  sec 

Mean  square 
difference 

0.60^  scc2 

0.44''scc2 

0.09*  scc2 

excluding  differences  la^er  than  4.0  sec 
*-  excluding  differences  larger  than  2.0  sec. 
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Table  4.  Cimiparisoo  of  residual  varianee  in  seconds^  for  selected  groups 
of  earthquakes,  for  both  sinde^vent  and  JHD/JED  relocations,  and  udng 
different  methods  for  calculating  travd  tiiMs.  The  residual  variance  is  the 
square  of  the  average  RMS  residual.  All  these  locations  utilized  input  data 
read  personally  by  the  author  or  his  close  colleagues. 


Region 

Single-Evem 

IASFEI91 

JHD/JED 

IASPEI91 

Single-Event 

J-B 

JHD/JED 

J-B 

Macquarie  Ridge 
(36  earthquakes) 

1.25 

0.83 

1.35 

0.71 

Bucaramanga 
with  pP  and  sP 
(26  earthquakes) 

5.02 

0.52 

• 

• 

Bucaramanga 
with  pP,  no  sP 
(26  earthquakes) 

4.97 

0.48 

• 

• 

Bucaramanga  -  P  only 
(26  earthquakes) 

4.20 

0.28 

1.30 

0.28 

Region 

Singte'Event 
"best"  model 

JHD/JED 
"^t"  model 

Single-Event 
layer  over 
h^space 

JHD/JED 
layer  over 
halfspace 

Efaie  -  P  and  S 
at  all  stations 
(21  earthquakes) 

0.23 

0.07 

0.23 

0.06 

Efate  -  P  and  S 

0.21 

0.04 

0.21 

0.04 

except?  only 
at  lai^  stations 
(21  earthquakes) 
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An  impGrtant,  fundamental,  and  presently  unanswered  question  is  what  might  cause  such  a 
penitent  souree  of  earthquake  activity  within  the  mande.  I^sumaldy  it  is  related  somehow  to  die 
subducdon  oi  the  Caribbean  or  Nazca  plates  beneath  South  America  (Pennington,  1981). 
Ifowever,  Schneider  et  al  (1987)  observed  that  P  and  T  axes  of  focal  mechanisms  determined  from 
the  temporary  local  netwoik  were  oriented  almost  randomly  ova-  the  focal  sphere  (Hgure  8).  This 
is  inconsistent  widi  the  usual  models  explaining  intermediate  focus  earthquake  mechanisms,  where 
they  are  due  to  failure  of  a  lithospheric  plate  due  to  bending,  <x  to  repeated  motion  along  a 
particular  weak  failure  surface. 

Nevertheless,  for  several  leastxis  the  Bucaramanga  nest  provides  a  convenient  source  for 
testing  and  evaluating  teleseismic  locadoo  algorithms.  Because  die  nest  is  so  localired  in  space,  for 
any  particular  data  set  or  location  method  the  nest  volume  ot  the  variability  in  event  locations  is  a 
measure  of  the  relative  accuracy  of  the  locatitxi  method  (Fiohlich  and  Kadinsky>C^de,  1993). 
Because  the  nest  occurs  well  beneath  the  crust,  teleseismic  P  phases  are  mme  likely  to  be  clear  and 
impulsive.  Thus,  the  variability  in  the  relative  locations  determined  is  attributable  largely  to 
systematic  errors  in  the  locatitxi  method,  rather  than  to  intrinsic  variability  in  the  locations  or  to 
erroneous  arrival  time  picks. 

Data  Analysis:  Bucaramanga's  locatitm  in  ncnrthem  South  America  allows  us  to  obtain 
good  station  coverage  in  the  northern,  southern,  and  eastern  azimuthal  sectors.  Using  the 
WWSSN  microfiche  seismogram  collection  at  (Columbia  University,  we  searched  to  find  P  arrivals 
at  21  stations  (Figure  9)  fOT  26  different  earthquakes.  Generally,  we  could  pick  P  arrivals  at  the 
majority  of  these  stations  for  earthquakes  having  body-wave  magnitude  mb  of  5. 1  or  greater. 

The  carefully  (fetermined  phase  arrivals  and  relocations  of  these  26  earthquakes  provided 
the  basic  teleseismic  data  set  for  this  study.  For  comparison  we  consider  two  other  teleseismically 
determined  set  of  locatimis;  locations  reported  directly  by  the  ISC,  and  relocations  performed  for 
tins  study  using  ISC-repoited  arrivals  at  a  restricted  set  of  stations  (Table  5  and  Figure  9).  We  also 
compare  these  locations  to  27  hypocenters  repmted  by  Schneider  et  al  (1987)  determined  using  a 
temporary  local  network. 

Because  the  nest  is  so  localized  in  space,  fOT  map  view  and  cross  sectional  comparisons  we 
present  pictures  of  the  minimum- volume  convex  polyhedrons  enclosing  the  hypocenters  to  display 
tire  relocations  (Figure  10).  Frohlich  (1992)  describes  in  detail  an  algorithm  for  determining  the 
minimum-volume  convex  bounding  polyhedron  for  a  group  of  hypocenters.  Essentially  the 
polyhedron  is  formed  by  finding  all  the  triangles  determined  by  triples  of  hypocenters  within  the 
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T,  B,  and  P  axes  for  Bucaramanga  focal  mechanisms. 


O  -  Harvard  CMT;  a  as  other  teleseismic;  +  =  Schneider  et  al.  (1987) 
T  axes  B  axes  P  axes 


"Sum”  Focal  Mechanisms 

fMl  ZMi/2:ilMilI 

Teleseismic  only  All  tnerhanijemg 


Hgure  8.  (T<^)  Summary  of  P,  B  and  T  axes  of  focal  mechanisms  for  Bucaramanga  earthquakes. 
Symbols  are:  octagons  -  Harvard  CMT;  triangles  -  other  teleseismically  determined  mechanisms 
reported  by  Pennington  (1981);  pluses  -  mechanisms  determined  from  temporary  local  network 
di^  by  Schneider  et  aL  (1987). 

(Bottom)  Cmqxrsite  focal  mechanism  determined  by  adding  together  moment  tensors  for  (bottom 
left)  Harvard  CMT  and  otlwr  teleseismic  eardrquakes;  and  (bottom  right),  all  available  mechanisms, 
each  normalized  by  dividing  by  the  scalar  sdsmic  moment 
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Table  5.  Hypocenters  in  the  Bucaramanga  nest  relocated  using  the  JHO 

method  and  P  PKP,  and  pP  times  arrival  times  read  personally  by 
the  author. 


Dale 

Hr  Min  Sec 

1 

Longitude 

®W 

Depth 

km 

rms 

sec 

Tmes 

23:36:14.27 

6.796 

72.921 

158.53 

0.79 

7/30/65 

7:20:10.15 

6.819 

72.935 

159.24 

0.49 

9/11/66 

17:38:  4.80 

6.831 

72.910 

157.47 

0.71 

3/21/67 

18:11:44.43 

6.776 

72.925 

159.91 

0.65 

7/29/67 

10:24:26.22 

6.832 

72.995 

166.56 

0.78 

5/7/68 

9:  0:29.98 

6.850 

72.960 

162.17 

0.68 

8/18/72 

20:58:19.17 

6.802 

72.997 

161.59 

0.74 

Amm 

16:26:10.71 

6.785 

72.993 

160.73 

0.89 

11 

4;  3:36.52 

6.802 

72.926 

160.70 

0.62 

9/5/74 

7:47:43.13 

6.840 

73.005 

168.33 

0.83 

3/13/76 

21:44:42.41 

6.806 

72.906 

161.06 

0.56 

4/1/76 

19:21:16.39 

6.839 

73.009 

158.06 

0.65 

3/23/77 

2:11:15.65 

6.801 

72.889 

161.38 

0.42 

7/2/78 

2:49:19.18 

6.743 

72.938 

164.40 

0.61 

1/14/79 

19:20:31.45 

6.742 

72.983 

161.70 

0.88 

3/11/79 

12:16:20.29 

6.805 

72.956 

161.28 

0.70 

8/15/80 

21:30:47.34 

6.749 

72.984 

160.42 

1.05 

12/13/81 

21:40:37.68 

6.788 

72.956 

157.49 

0.58 

5/19/83 

23:22:16.16 

6.825 

72.934 

161.05 

0.47 

8/29/83 

8:24:25.33 

6.806 

72.976 

156.77 

0.74 

11/4/84 

13:14:21.32 

6.804 

72.980 

161.60 

0.44 

6/26/85 

23:44:55.48 

6.799 

72.921 

159.74 

0.84 

12/3/85 

17:52:25.09 

6.777 

72.911 

159.85 

0.32 

3/31/86 

1:19:25.23 

6.846 

72.976 

156.06 

0.85 

6/29/86 

20:11:52.29 

6.859 

72.949 

165.05 

0.66 

10/31/86 

7:  5:54.05 

6.826 

73.001 

159.94 

0.35 
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Figure  10.  Boumling  polyhednms  for  four  sets  of  hypocentral  locations  in  the  Bucanunanga  nest. 
At  left  are  nu^  views  of  ^  upper  surfaces  of  polyhedrcxis;  at  right  are  cross  sectional  views  of  the 
fiont  side  of  the  polyhedrmis  as  viewed  from  die  east  The  four  sets  of  locations  are:  A)  Locations 
as  rqiorted  by  the  ISC  for  26  hypocenten  occuning  betwem  19&  and  1986;  B)  Relocations  of  the 
same  26  hypocenters  detenniii^  for  this  study  widi  a  weighted  least-squares  JHD  method  using 
phase  arrivals  reported  by  the  ISC  at  28  selected  stations;  C)  Relocations  for  the  same  26 
hypocenters,  relocated  with  a  weighted  least-squares  JHD  method  using  P,  PKP,  and  pP  data 
personally  reread  by  the  audior  at  18  selected  stations;  D)  Locations  for  27  hypocenters  reported  by 
Schneider  et  al.  (1987),  determined  from  a  temporary  local  network. 
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group,  then,  the  minimum  volume  crnivex  polyhedron  is  just  the  region  in  space  enclosed  by  these 
triangles. 

It  is  possible  to  ccmipaie  the  absolute  locations  for  the  26  earthquakes  which  possessed 
diiee  different  locations,  i.e.,  locations  reported  by  the  ISC,  lelocaticMis  determined  for  this  study 
using  phase  arrivals  reported  by  the  ISC,  and  relocations  determined  for  this  study  using  phase 
arrivals  read  by  the  author.  These  comparisons  (Figure  11)  indicate  that  the  median  discrepancy 
between  events  in  the  first  and  second  sets  of  locations  (reported  and  relocated  ISC)  was  D12  - 
17.3  km,  while  the  median  discrepancy  between  the  first  and  third  sets  (reported  ISC  and  this 
study)  was  D13  »  12.4  km. 

For  the  393  earthquake  arrivals  reread  especially  for  this  study,  38%  agreed  with  the  ISC 
reported  arrivals  within  0.5  sec  (Table  3).  However,  7%  differed  from  the  ISC  reported  arrivals 
by  1.0  seconds  or  mote.  The  mean  difference  between  phase  arrivals  read  by  us  and  repented  by 
the  ISC  was  0.39  seconds.  Moreover,  26%  of  the  phases  we  found  and  reread  were  not  repented 
to  the  ISC.  Also,  there  were  ISC-repe}rted  phase  arrivals  for  16%  of  the  phases  which  we  were 
unable  to  frnd  after  a  careful  search  oi  the  WWSSN  microfiche  seisme)grams. 

Using  the  phases  reread  specifically  for  this  study,  we  performed  relocations  of  the  26 
earthquakes  using  both  IASPEI91  and  Jeffreys-Bullen  (J-B)  travel  times,  and  using  both  single¬ 
event  and  7HD  location  algorithms  (Table  5).  Altiiough  the  variance  of  travel  time  residuals  (Table 
4)  obtained  using  IASPEI91  times  were  somewhat  larger  than  diose  obtained  using  J-B  times,  this 
was  mostly  because  we  did  not  have  J-B  tables  for  depth  phases  such  as  pP  and  sP.  Thus  these 
phases  did  not  constrain  the  relocations,  and  while  the  resulting  variances  were  small  the  nest  so 
obtained  was  far  less  localized  than  that  obtained  with  the  1ASPEI91  tables  (Table  6). 

The  volume  of  the  bounding  polyhedron  depended  strongly  on  exact  phases,  travel  time 
model,  and  method  of  location  used  (Table  5).  Bounding  polyhedrons  for  single-event  relocations 
had  volumes  of  30-200  times  larger  than  JHD  relocations.  This  was  caused  by  the  occurrence  of 
systematically  large  residuals  averaging  4-6  seconds  at  two  stations,  BHP  and  SJG,  which  also 
produced  variances  for  4-S  sec^  for  all  single-event  IASPEI91  relocations  (Table  4).  Thus,  for  the 
IASPEI91  relocations  using  JHD  reduced  the  variance  by  90%  or  more,  to  about  O.SO  sec^  and 
less.  The  effect  of  using  different  methods  to  calculate  travel  times  was  less  pronounced.  For 
example,  for  the  JHD  relocations  the  bounding  polyhedron  volume  was  2,135  km^  for  the  J-B 
relocations,  and  687  km^  for  the  IASPEI91  relocations  utilizing  P,  PKP  and  pP  phases. 
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Diff.  in  Latitude  (km)  Diff.  in  Latitude  (km) 


Diff.  in  Longitude  (km)  Diff.  in  Longitude  (km) 


Hgure  11.  Differences  in  locations  for  26  epicenters  located  by  three  different  methods;  1) 
locations  reported  by  the  ISC;  2)  locations  relocated  using  phases  reported  to  the  ISC;  and  3) 
locadtMis  determined  using  phases  reread  for  this  study.  Top  panels  summarize  differences 
between  the  first  two  sets  of  locations,  and  the  bottmn  two  panels  summarize  differences  between 
the  first  and  third  sets. 
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Table  6.  Comparison  of  volume  of  smallest  convex  polyhedron  enclosing 
relocated  groups  of  ^26  earthquakes  in  Bucaramanga,  Colombia,  and  21 
earthquakes  in  Efate,  Vanuatu,  for  both  single-event  and  JHD/JED 
relocations,  and  using  different  travel  time  tables.  All  these  locations 
utilized  input  data  read  personally  by  the  author  or  his  close  colleagues. 


Travel  time, 
phase  data 

Single-Event 

JHD 

Bucaramanga,  Colombia 

IASPEI91, 

P,  PKP,  pP  and  sP 

33,941  km3 

1,176  km3 

IASPEI91, 

P,  PKP  and  pP 

45,766  km3 

687  km3 

IASPEI91. 

P  and  PKP  only 

208,025  km3 

851  km3 

J-B,  P  only 

17,819  km3 

2,135  km3 

Efate,  Vanuatu 

"best"  model, 

P  and  S  at 
all  stations 

5,502  km3 

4,321  km3 

layer  ova* 
htdfspace,  P  and  S 
at  all  statitms 

7,227  km3 

5,022  km3 

"best"  model, 

P  and  S  except 

P  only  at  land 
stations 

8,55?  km3 

6,726  km3 

layer  over 
halfspace,  P  and  S 
except  P  only  at 
land  stations 

7,202  km3 

7,246  km3 
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C.  EfiM,Vanuatu 


Background:  Vanuatu  (formerly,  the  New  Hebrides)  is  a  nearly  linear  island  chain 
(Figure  12)  situated  at  the  convergent  boundary  between  the  Indo-Australian  and  Pacific  plates 
(Pascal  et  al.,  1978).  Near  the  island  of  Efate  the  convergence  rate  is  about  1 1  cm/yr  and  the 
direction  is  iqrproximately  normal  to  the  New  Hebrides  trench  (Pievot  et  al.,  1991).  About  TOO  km 
to  die  north  there  is  a  peculiar  region  of  shallow  bathymetry  known  as  the  d’Entrecasteaux  Ridge 
(DR)  which  clearly  affects  the  regional  tectonics  (C^ollot  and  Fisher,  1991;  Fisher  et  al.,  1991). 
The  ctdlision  d  die  DR  with  the  New  Hebrides  trench  appears  to  cause  (^atemary  uplift  and  tilting 
on  the  islands  of  Santo  and  Malekula  (Taylor  et  al.,  1980;  1987)  and  reactivation  of  normal  faults 
on  the  islands  of  Santo,  Maewo  and  Pentecost  (Isacks  et  al.,  1981).  Moreover,  while  the  New 
Hebrides  trench  is  clearly  evident  west  of  Efate  and  the  islands  to  the  south,  it  disappears  between 
Efate  and  Malekula,  and  the  western  parts  of  the  islands  of  Malekula  and  Santo  lie  along  the 
extenskm  of  the  New  Hebrides  Trench. 

The  DR  may  influence  the  ongoing  seismicity  of  the  central  Vanuatu  region  as  well.  Just 
west  of  Efate  near  n.S^’S,  168.0°E  there  is  a  highly  active  zone  of  shallow  seismic  activity  known 
as  the  Efate  nest  (Figure  12).  The  origin  of  the  Efate  nest  is  uncertain.  Chatelain  et  al.  (1986) 
suggested  this  area  forms  a  seismic  boundary  between  a  quiescent  asperity  to  the  north  and  a 
seismically  active  segment  in  the  south.  Prevot  et  al.  (1991)  observed  slow  seismic  velocities  for 
seismic  waves  traversing  die  Efate  nest  region.  They  suggest  that  the  Efate  nest  may  occur  where  a 
previously  subducted  remnant  of  the  DR  trails  to  the  south  of  its  bathymetric  expression  and 
interacts  with  the  overlying  plate. 

The  most  significant  uncertainty  in  the  interpretations  of  Chatelain  et  al.  (1986)  and 
Prevot  et  al.  (1991)  concerns  the  focal  depth  of  the  earthquakes  in  the  Efate  nest.  All  their  data 
comes  from  a  regional  seismograph  network  operated  between  1978  and  1989  at  up  to  23  sites  on 
various  islands  of  Vanuatu  by  the  Institut  Frangais  de  Recherche  Scientifique  pour  le 
Developpement  en  Cooperation  (ORSTOM).  Because  of  the  geometry  of  the  islands  forming 
Vanuatu  fix  the  ^metry  of  the  ORSTOM  land  mtwork,  it  is  difficult  to  obtain  reliable  depths  for 
earthquakes  occurring  northwest  of  Efate  and  south  of  Malekula. 

In  order  to  more  accurately  determine  earthquake  focal  depths,  between  June  and  August  of 
1989  ORSTOM  atxl  the  University  of  Texas  at  Austin  deployed  ocean  bottom  seismographs  (OBS) 
at  26  sites  within  this  region.  While  the  data  collected  by  the  OBS/land  network  was  not  originally 
collected  or  analyzed  for  the  present  study,  it  provides  a  unique  data  set  for  comparison  with  the 
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Figure  '  .  Map  the  central  New  Hebrides  island  arc,  with  major  islands  of  Vanuatu  labeled.  Box 
delineates  region  of  map  in  Figure  13  describing  1989  OBS/island  network.  Contours  are  in 
meters.  D'Entrecasteaux  Ridge  (labeled  "DR”)  lies  nenth  of  westward  bend  in  40(X)  m  contour 
occurring  west  of  Santo  and  Malekula.  Plus  symbds  (”+")  are  locations  reported  by  Roecker  et  al. 
(1988)  in  region  with  unusual  deep  focal  ^pths.  Circles  are  locations  reported  by  ISC  for 
earthquakes  with  arrivals  reported  by  OBS  array  in  Table  8.  Shaded  area  west  of  Efate  delineates 
the  Efate  nest-  a  regiem  of  intense  shallow  earthquake  activity  shidied  by  Chatelain  et  al.  (1986); 
within  the  shaded  area  between  1977  and  1988  the  ORSTOM  land  network  reported  30  or  mme 
eai^uakes  occurring  widiin  each  .02*^  x  .02'*  square,  with  each  location  constrained  by  at  least  7  P 
arrivals  and  4  S  arrivals. 
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two  teleseiinik  dam  sets.  It  also  provides  a  reasonable  test  ctf  the  ^;>i^catkMi  of  the  TexFlex 
locttioo  program  as  ^)|died  to  local  network  data. 

Data  Anafysis:  The  field  esqteriment  consisted  of  two  deployments  of  13  digital  OBS 
utttnnnents  each  (Hgme  13)  between  22  June  and  19  August,  1989,  for  a  total  dei^oyment  time  of 
42  days.  The  OBS  stations  in  this  study  operated  in  a  triggered  mode  as  described  by  Irohlich  et 
al.  (1990),  widi  the  trigger  dqtending  on  dm  ratio  oi  absolute  signal  amplitudes  determined  over  a 
short-term  (1.6  seccMxl)  and  a  long-term  (411  second)  interval.  A  continuously  operating  circular 
memory  buffer  held  6- 10  seconds  of  data,  so  diat  wimn  the  OBS  triggered  a  t^  drive  recorded  dm 
signal  in  the  buffer  as  well  as  the  signal  fdlowing  the  trigger.  All  OBS  had  a  vertical-component 
4.5  Hz  seismometer.  In  additirm,  all  but  two  instruments  had  two  4.5  Hz  hraizraital-component 
seismoinetBrs.  One  remaining  instrument  had  a  hydrophone  and  one  4  J  Hz  horizontal-ctHnpramnt 
seismometer  (OB5  and  021);  dm  last  instrument  had  a  2.0  Hz  vertical-componrot  seismometer  and 
a  10  Hz  horizontal  comptment  seismometer  (OBS  and  026).  The  sampling  interval  for  recording 
was  25.056  ms,  widi  data  low  pass  filtered  with  a  cut-off  frequency  of  15  Hz  prior  to  recording. 

The  data  also  included  records  fmn  niim  land  stadrais  in  the  ORSTOM  regicmal  network 
(Figure  13).  These  stations  all  have  vertical  component,  4.5P  Hz  seismometers.  Data  are 
telemetered  to  a  base  staticm  where,  if  enough  stations  trigger,  data  are  recorded  cm  magnetic  uqm 
and  by  a  multichannel  oscillograph  system.  Experience  suggests  that  the  system  records  all 
earthquakes  in  dm  study  area  with  magnitudes  Ml  larger  than  about  2.7  (Chatelain  et  al.,  1986). 

Altogether  the  land  and  OBS  networks  recorded  about  40,000  blocks  of  data  during  this 
ejqmriment  For  dm  initial  data  processing  we  examined  each  block  individually  to  identify  seismic 
signals.  For  the  events  identified  as  earthquakes  occurring  within  the  Efate  nest  (Figure  14),  two 
different  sdsmologists  independently  picked  the  times  of  P  and  S  phases,  and  evaluated  the  quality 
of  the  arrival  time  picks  by  grading  them  from  1  to  4,  corresponding  approximately  to  uncertainties 
of  <0.1  sec,  <0.3  sec,  <1.0  sec,  and  >1.0  sec,  respectively.  For  5,991  arrivals  assigned  pick 
qualities  1  and  2,  the  mean  difference  in  picked  times  was  0.21  sec  (Table  3),  and  90%  of  all 
differences  were  less  than  0.55  seconds. 

Of  3443  seismic  events  recorded  during  the  field  study,  28  were  airgun  shots  that  we  fired 
during  the  field  program,  18  were  regional  or  teleseismic  earthquakes  reported  by  the  International 
Seismological  Centre  (Figure  12),  and  74  were  observed  only  at  stations  in  the  land  network 
without  triggering  any  OBS.  We  further  classified  2563  of  these  events  as  unlocatable  because 
they  possessed  P  phases  of  pick  quality  3  or  better  at  2  or  fewer  stations,  or  at  3  stations  with  no  S 
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Hgure  14.  Histogram  of  difference  in  P  and  S  phase  anival  times  as  indeed  indq)midently  by  two  individuals. 


l^uttes  of  pidc  quality  3  or  better.  Unfortunately,  none  of  the  earthquakes  reported  by  the  ISC 
were  well-iecotded  enough  by  the  OBS  network  to  locate  independently. 

For  the  remaining  760  earthquakes,  we  obtained  preliminary  locatkms  using  TexFlex  in  the 
single  event  model  (Figure  15).  Our  "best"  velocity  model  was  similar  to  those  determined  by 
Prevotetal.  (1991)  for  the  region  beneath  and  just  west  of  Efate.  For  arrivals  at  the  OBS  stations, 
the  "best”  model  had  three  layers  of  thickness  10  km  overlying  a  halfspace,  with  velocities  of  S.O 
km/sec,  7.1  km/sec,  7.3  km/sec,  and  8.0  km/sec.  At  the  land  stations,  the  layer  thicknesses  were 
10  km,  20  km,  and  10  km,  and  the  velocities  were  5.2  km/sec,  6.8  km/sec,  7.6  km/sec,  and  8.0 
km/sec.  For  comparison,  we  also  poformed  all  relocations  assuming  a  simpler  model  consisting 
of  a  layer  of  velocity  S.S  km/sec  overlying  a  halfspace  with  velocity  8.0  km/sec.  For  both  the 
"best"  aixl  the  "layer  over  a  halfspace"  models  we  assumed  that  Vp/Vs  was  1.78.  The  remaind^  of 
our  analysis  ctmcems  evaluation  of  the  subset  of  these  760  earthquakes  which  occurred  within  the 
Efate  nest 

To  obtain  ^ecise  relative  locations  for  events  within  and  near  the  Efate  nest  we  first  appired 
die  JHD  method  to  an  especially  well-recorded  subset  of  21  earthquakes  in  ortter  to  determine 
stadcHi  corrections.  At  least  four  OBS  stations  and  two  land  stations  recorded  each  of  these 
earthquakes,  with  an  azimuthal  gap  of  230°  or  less.  Smne  experimentation  revealed  that  the  JHD 
station  corrections  depended  strongly  on  the  particular  combinations  of  land  stations  data  used. 

We  finally  chose  to  relocate  using  both  P  and  S  data  from  OBS  stations  while  using  P 
arrival  data  only  from  the  land  stations  DVP,  RTV,  MBV,  and  NGA.  With  this  scheme  the  (S  -  P) 
times  at  three-component  OBS  stations  constrained  focal  depths  of  events,  while  P  times  at  the 
more  distant  vertical-component  land  stations  influenced  the  lateral  position  of  the  events.  The 
resulting  station  corrections  were  positive  at  the  land  stations  in  the  east  (NGA,  MBV  and  RTV), 
quite  negative  at  four  stations  near  the  center  of  the  network  (OB3,  OB4,  OB7  and  025),  and  quite 
positive  for  two  of  the  OBS  stations  at  the  western  edge  of  the  netwotic  (OlO  and  027).  Applying 
these  station  conections  reduced  the  variance  for  the  21  especially  well-tec(xded  earthquakes  from 
0.21  sec^  to  0.04  sec^  (Table  4). 

While  the  application  of  station  conections  did  reduce  the  variance  of  travel  time  residuals, 
it  had  only  a  minor  effect  on  the  localization  of  the  group  of  21  events,  as  measured  by  the  volume 
of  the  bounding  polyhedron  (Table  6).  Also,  there  was  no  systematic  difference  between  using  the 
"best"  velocity  model  and  the  layo'-over-halfspace  model. 

Hiuilly,  to  evaluate  the  overall  spatial  gecxnetry  of  the  events  recortted  from  within  the  Efate 
nest,  we  applied  the  station  corrections  obtained  above  to  relocate  a  larger  set  of  92  earthquakes 
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Hgure  15.  Map  of  preliminary  locations  of  earthquakes  recorded  by  OBS  stations  in  this  study. 
Otises  ("X")  and  pluses  ("+")  are  events  recorded  in  the  first  and  second  OBS  deployments, 
respectively.  Map  tmly  includes  locations  obtained  using  at  least  4  P  and  1  S  observation  and 
having  an  azimuthal  gap  of  330P  or  less.  The  box  delineates  the  region  mai^)ed  in  Figure  16. 
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occurring  near  Efate  (Figure  16).  Each  of  locadcMis  utilized  at  least  4P  and  2S  arrivals  and 
possessed  an  azimuthal  giq>  smaller  than  27(P. 

3.  EVALUATION  OF  SOURCES  OF  LOCATION  UNCERTAINTY 
A.  Background 

In  this  course  of  the  study  we  determined  earthquake  locations  using  a  broad  variety  of 
fairly  standard  methods,  i.  e.,  for  travel  times,  we  utilized  the  IASPEI91  and  J-B  models  as  well  as 
times  detetmined  from  tracing  rays  in  flat-earth  models;  for  station  cmtections  we  employed  both 
single-event  relocation  methods  as  well  as  JHD  and  JED  methods.  The  events  themselves  included 
shallow  eartlKiuakes  recorded  only  by  teleseismic  stations,  shaUow  earthquakes  recorded  tmly  by 
local  stations,  and  intermediate  focus  earthquakes  recorded  by  both  regional  and  teleseismic 
statkms.  For  this  study  the  flexible  character  of  the  TexFlex  relocation  {uogram  (Frohlich,  1993) 
was  essential,  as  it  allowed  us  to  carry  out  most  of  the  different  relocations  above  simply  by 
changing  a  single  parameter  in  the  data  input  flies. 

The  objective  in  this  section  is  u>  compare  and  evaluate  the  locations  in  different  geography 
regions  to  reach  some  overall  conclusions  about  the  sources  of  location  uncertainties  for  events 
located  using  travel-time  based  methods.  Fot  this  study  we  have  three  different  ways  to  evaluate 
die  quality  of  locations;  tiiese  are: 

1)  Dependence  of  the  variance  of  travel  time  residuals  on  the  method  of  location  and  the 
selection  of  data  (Table  4); 

2)  Difference  in  geographic  locations  as  determined  by  independent  agencies  (Figures  6  and 
11);  and 

3)  Dependence  of  volume  of  the  bounding  polyhedron  oi  the  method  of  location  and  selection 
of  data  (Table  6  and  Figure  10). 

For  the  purposes  of  discussion  of  travel  time  variances,  we  shall  assume  that  location 
uncertainties  contribute  approximately  additively  to  the  variance  Vtot  of  travel  time  residuals  for 
groups  of  locations  determined  by  a  single-event  method.  It  is  reasonable  to  assume  that  errors 
contribute  additively  as  the  square  of  residuals  since  all  the  location  algorithms  used  are  least- 
squares  algorithms,  i.e.,  they  minimize  the  sum  of  squares  of  residuals. 

Thus,  to  formulate  a  crude  model  of  the  location  uncertainty  we  posit  that: 

Vtot  *  Vdata  +  Vstations  +  Vvetocity-model  +  Vothcr  =  (RMS)^ 
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Figure  16.  Miqp  of  relocations  for  earthquakes  near  Efate.  Crosses  and  pluses  are  as  in  Hgure  IS; 
heavy  daik  symbols  are  for  21  especially  well-recorded  earthquakes  used  to  (tetermine  station 
corrections;  lighter  symbols  are  71  moderately  well-recorded  events  (see  text).  Lines  maticed  "XL 
-  XR"  and  "YL  -  YR"  show  locations  of  cross  sections  in  Hgure  17. 
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and  estiinaie  the  relative  sizes  of  each  of  the  components  of  Vtot  fw  the  different  location  regimes. 
Here,  Vdbtu  represents  variance  attributable  to  incorrect  identification  of  phase  arrivals  or 
misreading  phase  arrivals;  Vgiatioos  represents  variance  attributable  to  near-station  structure,  ch'  to 
path  effects  which  strongly  influence  a  single  station;  V velocity-model  represents  variance  which 
affects  all  stations,  and  is  attributable  to  choosing  an  inai^ropriate  model  for  relocating  the  group 
of  events;  and  Vother  represents  other  sources  of  variance,  or  sources  of  variance  that  should  be 
assigned  to  Vdam.  ^sia&oas*  or  Vvelocity-inodd  but  which  our  analysis  is  unable  to  detect  While  this 
model  clearly  does  not  account  properly  feu*  all  sources  of  error  (e.g.,  it  wrongly  suggests  that 
earthquake  locations  with  zero  residuals  possess  no  locatitm  error),  it  is  a  useful  construct  for 
comparing  uncertainties  of  groups  of  earthquake  locations  determined  by  a  variety  of  different 
least-squares  methods. 

B .  Misreading  or  Incorrect  Identification  of  Phase  Arrivals 

One  approach  to  estimating  the  variance  Vdata  due  to  misreading  or  mispicking  of  phase 
data  is  to  evaluate  the  differences  in  phase  arrival  times  as  picked  independently  by  different 
readers.  If  we  imagine  a  hypothetical  situation  where  both  the  velocity  model  used  to  determine 
travel  times  and  the  station  corrections  were  known  exactly,  then,  the  only  source  of  variance 
would  be  misidentification  of  phases  or  picking  errors.  If  two  readers  each  pick  times  so  that  their 
pick  errors  have  variance  Vi  and  V2  respectively,  tten  the  difference  in  picks  will  have  variance  V 1 
+  V2.  Thus,  if  Dmean-squate  is  the  mean  square  difference  in  travel  time  picks  for  two  independent 
readers,  and  if  Vi  =  V2,  an  approximation  for  the  true  variance  would  be: 

Vdata = Dmean-squarc/  2 

Since  similarly  trained  readers  evaluating  the  same  data  are  likely  to  make  the  same  errors  in 
reading  or  phase  identification,  half  the  mean  square  difference  in  picks  (Table  3)  probably 
underestimates  the  true  size  of  Vdata- 


For  Macquarie,  Dmean-square  is  0.60  sec^,  suggesting  that  Vdata  is  at  least  0.30  sec^  or 
greater.  Fot  Bucaramanga,  Dmean-square  is  0.44  sec^  suggesting  that  Vdata  is  at  least  0.22  sec^  or 
greater.  For  Efate,  Dmean-square  is  0.09  sec^,  suggesting  that  Vdata  is  0.045  sec^  or  larger. 
Represented  as  fractions  of  Vtot  as  estimated  from  the  single-event  variance,  these  estimates  of 
Vdata  represents  24%  of  Vtot  in  Macquarie,  4%  in  Bucaramanga,  and  20%  in  Efate  (Table  7). 
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C.  Nw-^iuhH  or  Station-Path  ^cts 


Since  the  JED  and  JHD  methods  were  specifically  designed  to  adjust  event  relocations  for 
stadon-assignaUe  scNirces  of  location  uncertainty,  a  good  estimate  of  the  size  of  Vcuttions  is  the 
difference  in  variance  for  event  groups  relocared  with  single-event  and  with  JED  or  JHD  methods. 
For  Macquarie  Ridge  die  difference  (Triile  4)  in  variance  between  single-event  and  JED  relocations 
was  0.42  sec^  for  the  IASPEI91  model,  or  33%  of  the  total  variance  of  1.25  sec^  (Table  7).  For 
the  J-B  model  the  difference  was  0.64  sec^,  or  47%  of  the  total  variance  of  1.35  sec^.  For 
Bucaramanga,  the  IASI^91  variance  difference  was  about  4.50  sec^,  m'90%  of  the  total  variance 
ifriien  the  relocations  were  performed  both  with  and  without  sP  phases.  Nearly  all  this  variance 
difference  was  attributable  to  large  systematic  residuals  at  two  stations.  SJG  and  BHP.  For  Efate, 
Vanuatu,  relocations  using  P  and  S  phases  at  all  stations,  the  variance  difference  was  0. 17  sec^,  or 
73%  of  the  total  single-event  variance  of  0.23  sec^. 


Table  7.  Approximate  percentage  of  travel  time  residual  variance  attributable  to 
various  sources  in  each  of  the  three  geographic  regions  studied.  Sources  are: 
Vdaia  ’  variance  due  to  misreading  or  mispicking  phase  arrival  times;  Vgtatioa-path 
•  variance  due  to  near-station  structure  or  path  effects  unique  to  a  Articular 
station;  Vvciocity-nodei  *  variance  attributable  to  incorrect  choice  of  velocity 
model;  Vother  *  estimate  of  remaining  variance.  Sum  of  estimates  do  not 
necessarily  add  up  to  100%  because  all  estimates  are  highly  approximate. 


Regkm 

Vdna 

Vstatioo-paih 

^velocity-model 

Voiher 

Macquarie  Ridge 

24% 

30-50% 

10% 

16-36% 

Bucaramanga 

4% 

90% 

- 

6% 

Efate,  Vanuatu 

20% 

73% 

10% 

<0% 

Thus  we  conclude  that  the  station-dependent  contribution  to  location  variance  is  of 
considmble  impmtance.  In  situations  such  as  Macquarie  and  Efate,  Vanuatu  where  the  degree  of 
velocity  inhtnnogeneity  is  apparently  more-or-less  "normal"  it  is  still  possible  to  reduce  variance  by 
30-75%  simply  by  detmmining  statitm  corrections.  In  situations  such  as  Bucaramanga  where  there 
may  be  considerable  lateral  heterogeneity  due  to  the  presence  of  subducted  lithosphere,  the  variance 
leductimi  may  be  90%  or  more. 
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D.  Velocity44odel  Dependent  Effects 


To  estimate  the  velocity-model-dependent  souices  of  location  variance  Vveiocity-model 
compare  variances  of  event  groups  relocated  the  JED  ot  JHD  method  using  diffoent  velocity 
models.  For  Macquarie  Ridge,  the  variance  difference  (Table  4)  between  IASPEI91  aiul  J-B 
relocations  was  0.12  sec^  fm*  JED  relocations,  corresponding  to  about  10%  of  the  total  single¬ 
event  variance  of  1.25  sec^  (Table  7).  For  Bucaramanga,  there  was  no  variance  difference  for  P- 
<Mily  J-B  relocations  using  the  J-B  or  IASPEI91  mottels.  For  Efate,  Vanuatu,  the  variance 
difference  between  our  preferred  model  and  the  simple  layer-over-halfspace  model  was  only  0.02 
sec^,  or  about  10%  of  the  total  single-event  variance  of  0.23  sec^. 

Thus,  in  all  three  regions  the  velocity-model  dependent  contribution  to  location  variance  of 
about  10%  or  less  was  smaller  than  either  the  station-dependent  or  tiie  data-dependent  cmtribution. 
However,  we  note  that  in  the  present  study  we  were  careful  in  all  three  geographic  areas  to  evaluate 
(mly  events  for  which  we  had  reasonable  azimuthal  control.  For  events  with  poor  azimuthal 
control  the  velocity-model-dependent  location  errors  may  be  "hidden,"  since  the  least-squares 
method  can  then  just  adjust  all  the  locations  to  accommodate  differences  in  velocity  model. 
Qeariy,  when  azimuthal  control  is  poor,  a  proper  velocity  model  is  essential  for  obtaining  good 
locatitHis. 

£.  Other  Effects  -  Discussion 

In  this  study  Vother  is  the  variance  remaining  after  we  have  estimated  Vdaia,  Vstati(ms>  and 
Vveiocity-model-  For  Macquarie  Ridge,  our  approximate  estimates  above  (Table  7)  suggest  that 
Vother  is  16-36%  of  the  Vtot,  i.e,,  that  either  16-36%  of  the  variance  is  attributable  to  other 
sources,  or  that  we  have  underestimated  the  importance  of  one  or  more  of  the  sources  of  variance. 
Fen:  both  Bucaramanga  and  Efate,  the  estimates  of  Vother  ttre  so  small  that  they  are  indistinguishable 
from  zero. 

Incidentally,  the  size  of  differences  in  geographic  location  (Figures  6  and  11)  for  event 
locations  reported  by  independent  sources  is  roughly  consistent  with  size  of  travel  time  variances. 
In  particular,  in  both  Macquarie  and  Bucaramanga  the  geographic  differences  were  typically  about 
15  km,  whereas  the  variances  for  single-event  J-B  locations  were  1.35  sec^  and  1.30  sec^, 
respectively.  A  variance  of  1.3  sec  corresponds  to  an  RMS  travel  time  residual  of  1.15  sec.  Now, 
dT/dA  for  phases  between  30®  and  90®  ranges  from  8.7  sec/degree  to  4.6  sec/degree,  or  between 
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.08  Mc/km  and  .04  sec/km.  Thus,  geognq>luc  differences  15  km  cmrespond  roughly  to  travel 
time  residuak  (rfO.6  sec  to  1.2  sec. 

However,  the  data  concerning  the  effect  of  location  method  on  the  volume  of  bounding 
polyhedrons  demonstrates  that  travel  time  variance  itself  gives  an  incomplete  picture  oi  location 
uncertainty.  The  Bircaramanga  relocations,  possessing  reasonably  good  distribution  of  stations 
and  phases  which  are  generally  impulsive,  represent  about  as  good  accuracy  as  one  can  t^tain  by 
reading  analog  seismograms  at  teleseismic  stations.  Yet,  in  Bucaramanga  the  JHD  relocations 
using  P-phases  only  had  a  significantly  smaller  variance  than  those  with  P-  and  pP  phases  (0.28^ 
sec  compared  to  0.48^  sec,  see  Table  4,  however,  the  bounding  polyhedron  volume  was  slightly 
analler  for  ibt  relocations  with  P-  and  pP  (851  km^  vs  687  km^  see  Table  5).  Apparently  the  use 
of  dq)th  phases  better  constrains  the  hypocenteis  even  thou^  it  increases  the  variance. 

Also,  these  exercises  demonstrate  tiiat  adding  more  data  can  sometimes  actually  degrade  the 
quality  of  locations.  For  example.,  in  Bucaramanga  the  relocations  that  incmporated  sP  phases  had 
both  larger  bounding  polyhedrons  (1,176  km^  vs  687  km^)  and  larger  variances  (0.52  sec^  vs 
0.48  sec^)  than  did  locations  using  only  P,  PKP  and  pP. 

Finally,  the  comparison  of  phase  arrival  picks  read  specifically  for  this  study  with  those 
reported  by  the  ISC  (Table  3)  indicated  that  when  a  strong  phase  existed  and  was  reported,  the  time 
picks  were  in  remarkably  good  agreement,  usually  within  less  than  a  half  second.  A  more 
troubling  result  is  that  in  a  substantial  fraction  of  cases  either  no  phase  was  reported  (15%  and  26% 
for  Macquarie  and  Bucaramanga),  or  a  phase  was  reported  at  a  time  when  we  were  unable  to 
identify  any  significant  arrival  on  the  microfiche  (24%  and  16%  fOT  Macquarie  and  Bucaramanga). 
This  suggests  that  blunders  are  a  significant  source  of  error  in  relocations  using  ISC  reputed 
{tiiases.  This  will  have  only  minor  effects  on  research  projects  where  there  is  ample  data  and  the 
results  effectively  average  all  readings,  such  as  determining  the  hypocenters  cf  large,  well-recorded 
earthquakes.  However,  these  blunders  can  adversely  affect  research  projects  where  there  are  few 
data,  such  as  the  relocation  of  small,  poorly  recorded  earthquakes  or  explosions.  Similarly,  they 
could  be  crucial  for  projects  where  the  results  depend  strongly  on  a  few  anomalous  individual 
phase  readings,  such  as  travel  time  tomography  in  regions  where  there  is  significant  lateral 
heterc^eneity. 

One  implication  of  this  study  is  that  there  are  inherent  weaknesses  in  any  location  method 
which  relies  solely  on  travel  times.  In  particular,  in  specific  cases  we  have  demonstrated  that 
whenever  data  is  sparse  large  location  errors  are  likely  to  arise  from  the  effects  of  misread  or 
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unread  phase  enivals.  Only  when  the  number  of  observing  stations  is  large  will  misread  phases  be 
identified  as  "outliers"  in  the  pattern  oi  residuals. 


One  possible  pnnnising  apjnoach  that  works  even  when  the  number  of  observing  stations 
is  small  is  to  combine  travel  time  analysis  with  wavefram  modeling.  While  excellent  methods  have 
been  available  for  some  time  tur  modeling  long  period  seismograms  at  teleseismic  distaiKxs,  there 
have  recently  been  significant  advances  in  the  fcnward  modeling  of  synthetic  seismograms  at 
legicHial,  broadband  stations  (e.g.,  Zhao  and  Helmberger,  1991a;  1991b;  1993a;  1993b).  Almost 
invariably  forward  modeling  will  constrain  source  depth  more  effectively  than  travel  time  analysis. 
Moreover,  it  is  possible  to  locate  an  event  if  there  is  even  one  well-recorded  digital,  three- 
component  record.  While  this  is  not  always  available,  forward  modeling  to  match  available 
seismograms  should  place  enough  constraints  on  the  locations  so  that  one  can  identify  blunders  in 
]rfiase  arrival  identification. 

F.  Summary  and  Conclusions 

1.  We  have  utilized  a  specially  developed  relocation  program  -  TexFlex  -  to  evaluate  three 
groups  of  seismic  events  representing  markedly  different  location  scenarios;  geographically  these 
events  are  situated  along  the  Macquarie  Ridge;  in  the  mantle  beneath  Bucaramanga,  Colombia;  and 
east  of  Efate,  Vanuatu. 

2.  We  here  develop  an  approximate  model  to  evaluate  seismic  location  uncertainties  by 
partitioning  variance  in  travel  time  residuals  between  three  sources:  variance  caused  by  errors  in 
identifying  and  picking  phases;  variance  which  can  be  reduced  by  assigning  station  corrections; 
and  variance  produced  by  using  an  inccnnplete  or  inapprqrriate  velocity  model  for  relocation. 

3.  Fot  events  in  all  three  geographic  regions  the  largest  identifiable  source  of  variance  is 
that  associated  with  station  corrections,  representing  30-90%  of  the  total  variance  for  single-event 
relocations  with  no  station  correction. 

4.  The  next  largest  source  of  variance  was  that  associated  with  phase  misidentification  or 
mispicking,  representing  4-24%  of  the  total  variance.  For  well-recorded  phases  there  was 
remarkably  little  disagreement  in  the  phase  picks  as  determined  by  different  readers,  as  the 
differences  in  arrival  times  averaged  only  0.2-0.4  sec.  A  more  troubling  result  was  that  for  the 
teleseismic  data,  there  were  often  phase  arrivals  reported  at  times  where  we  could  identify  no 
arrival,  tx’  clear  arrivals  for  which  no  phase  was  reported.  This  suggests  that  in  practice  picking 
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blunders,  grossly  misidendfied  phases,  and  unreported  phases  are  mwe  significant  sources  of 
error  for  teleseismic  locations  than  is  misreading  of  times  for  strong  phases.  Automatic  picking  of 
phase  arrivals  might  help  avoid  these  problems,  although  this  may  be  unduly  optimistic. 

5.  In  this  study  the  smallest  source  of  variance  was  associated  with  an  inappropriate 
velocity  model.  This  source  of  variance  was  almost  negligible  for  two  of  the  three  geographic 
regions,  which  were  groups  of  events  relocated  using  station  cmrections  and  with  data  having 
good  azimuthal  coverage.  This  suggests  that  if  the  objective  is  to  determine  good  relative 
locations,  when  the  azimudial  distribution  of  stations  is  good,  it  is  often  unnecessary  to  determine 
very  detailed  velocity  models.  If  azimuthal  distribution  of  stations  is  poor,  accurate  locations 
require  a  proper  velocity  model. 

6.  A  previously  developed  way  of  evaluating  location  methods  is  to  determine  die  volume 
of  the  smallest  bounding  polyhedron  of  groups  of  events  known  to  originate  from  a  localized 
source  region  -  the  assumption  is  then  that  smaller  polyhedron  volumes  are  associated  with  more 
accurate  methods  for  determining  relative  locations.  This  approach  reveals  some  inadequacies  in 
die  residual  variance  analysis. 

7.  Reducing  residual  variance  is  not  always  associated  with  smaller  polyhedron  volume. 
Mtneover,  adding  more  data  sometimes  actually  degrades  the  locations,  presumably  when  the 
additional  data  is  of  relatively  poor  quality.  This  implies  that  when  data  is  sparse  the  qitimal 
location  strategy  for  a  particular  region  should  be  to  obtain  data  firom  a  relatively  few  good  stations 
(say,  6  stations),  selected  for  high  quality  and  for  azimuthal  coverage. 


5.  RESULTS  OF  REGIONAL/TECTONIC  INTEREST 
A.  Macquarie  Ridge 

Two  features  of  earthquakes  in  the  Macquarie  Ridge  region  are  worthy  of  special  note. 
First,  the  epicenters  (Figure  3)  seem  to  occur  in  several  distinct  clusters,  possibly  associated  with 
the  occurrence  of  a  large  earthquake.  These  clusters  are  separated  by  regions  where  activity  is  low 
or  absent  In  particular,  these  include: 

1)  A  clusto'  occurring  between  about  46.5°S  and  48*’S,  approximately  in  the  region  of  the 
me  =  7.3  thrust  earthquake  of  12  October,  1979; 
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2)  A  cluster  occioring  betwe«)  about  49**S  and  50°S,  including  the  aftershodc  regitm  the 
Mw  *  7.6  strike>slip  earthquake  of  25  May.  1981; 

3)  One  or  more  clusters  occurring  between  about  51**S  and  SS°S,  including  the  aftershock 
legion  of  the  great  Mw  ^  8.1  strike-slip  earthquake  d  23  May,  1989; 

4)  A  cluster  occurring  between  about  S8.S**S  and  60°S,  near  the  epicenter  of  the  Mw  »  7.4 
earthquake  of  3  September,  1987. 

An  important  question  concerns  whether  the  regions  of  low  activity  between  the  major 
clusters  represent  perennial  regions  of  low  activity;  or  whether  instead  they  represent  the  sites  of 
future  great  earthquakes.  Although  the  locations  of  eardiquakes  prior  to  1964  (Figure  1)  are  poorly 
determined,  none  with  magnitudes  larger  than  7.2  (x  so  have  reported  locations  falling  in  the  gaps 
between  the  clusters  defined  by  the  post-1964  data. 

Anothm*  notable  feature  of  Macquarie  seismicity  south  (tf  47^*8  is  there  is  little  variability  in 
focal  mechanism  type  if  one  ignores  principal  axes  that  are  poorly  determined  because  they  lie 
along  the  equatorial  plane  of  earthquakes  with  large  non-double-couple  components.  That  is,  the 
majority  of  well-determined  mechanism  axes  (Figure  2)  are  nearly  horizontal  P  axes,  or  nearly 
vertical  B  axes,  as  expected  for  strike-slip  earthquakes.  Thus  the  presence  of  a  number  of  (}MT 
with  significant  non-double-couple  components  seems  to  account  for  much  of  the  alleged 
variability  in  mechanism  type  that  had  been  previously  noted  by  Ftohlich  et  al.  (1989). 

An  important  unresolved  question  concerns  the  cause  of  the  relatively  lar^  numbo*  of  non¬ 
double-couple  earthquakes  observed  in  this  region.  In  particular,  are  the  non-double-couple 
compcments  spurious,  coming  about  merely  because  the  region  is  so  remote  frmn  most  broadband 
digital  seismic  network  stations?  Or  are  the  non-double-couple  mechanisms  real,  related  in  some 
way  to  the  complexity  tectonics  of  this  region,  perhaps  caused  by  the  fact  that  the  pole  of  rotation 
for  Pacific-Australian  motion  is  so  close  to  the  Macquarie  Ridge? 

B .  Bucaramanga,  Colombia 

The  present  study  is  by  no  means  the  first  to  demonstrate  the  remarkable  character  of  the 
earthquake  activity  comprising  the  Bucaramanga  nest  Our  relocations  do  find  the  dimension  (~1S- 
20  km)  and  volume  (~700  km^)  of  the  nest  to  be  smaller  than  any  previous  study  depending  solely 
on  teleseismic  data.  Presumably  this  is  becau%  of  the  care  we  have  taken  in  using  only  data  read 
by  the  author,  and  because  we  have  relocated  using  a  set  of  stations  selected  K>  give  good  azimuthal 
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covenife  of  the  scnnce  region.  However,  these  results  are  still  considerably  less  localized  than 
those  found  by  Schneider  et  al.  (1987)  using  data  fiom  a  tempocary  local  network. 

One  important  difference  between  our  results  and  those  of  Schneider  et  al.  (1987)  concerns 
die  distribution  of  focal  mechanisms  for  earthquakes  within  the  Bucaramanga  nest  (Figure  8). 
When  we  consider  only  mechanisms  determined  from  teleseismic  observations,  we  find  that  the 
majority  of  the  mechanisms  possess  T  axes  which  dip  eastward  at  about  45^,  P  axes  which  dip 
westward  at  about  45**,  and  B  axes  which  are  approximately  north>soudi  and  horizontal.  Thus,  the 
resulting  "sum**  mechanism  has  one  approximately  htxizontal  nodal  plane,  ami  (me  ap[noximaiely 
vertical  mxlal  plane  striking  nordi-south. 

In  ccmtrast,  Schneider  et  al.'s  (1987)  mechanisms  determined  from  l(x:al  netw(xi(  data 
possess  principal  axes  which  are  scattered  over  almost  the  entire  focal  sphere.  We  can  only 
(xmcluite  that  either  there  is  a  genuine  difference  between  the  mechanisms  of  large  and  small 
earthquakes  within  the  nest,  or  that  there  was  scmie  unknown  error  in  the  data  or  methcxl  used  to 
determine  the  mechanisms  of  the  smaller  events.  In  any  case,  if  there  was  an  error  it  may  no 
longer  be  necessary  to  search  for  a  physical  explanation  of  the  Bucaramanga  nest  which  the 
extreme  range  of  focal  mechanisms  noted  by  Schneido'et  al.  (1987). 

C.  Efate,  Vanuatu 

Figure  16  and  the  accompanying  cross  sections  (Figure  17)  clearly  demonstrate  several 
important  characteristics  of  the  Efate  nest.  First,  all  but  one  of  the  hypocenters  has  a  focal  depth 
exceeding  2S  km.  Activity  is  absent  at  ncnmal  crustal  depths,  but  quite  intense  between  about  30 
km  and  45  km.  The  network  geometry  is  not  ideal  for  l(x:ating  deeper  events  (e.g.,  see  McLaren 
and  Frohlich,  1985),  but  those  visible  presumably  occur  within  an  eastward-dipping  Wadati- 
Benioff  zxmc. 

Sec(md,  the  hyp(x:enters  making  up  the  region  of  most  intense  activity  between  30  km  and 
50  km  depth  form  a  diffuse  cloud,  which  at  first  glance  contains  no  obvious  structure.  However, 
more  careful  scrutiny  suggests  that  this  cloud  includes  a  westward-dipping  planar  structure. 
Indeed,  50  of  the  92  hyp(x»nters  in  Figures  16  and  17  lie  within  3  km  of  a  planar  surface  having  a 
westward  dip  of  33**  and  strike  12**  W  of  N,  delineated  by  arrows  in  Figure  18. 

However,  this  planar  surface  is  clearly  not  parallel  to  either  ncxlal  plane  (Figure  18)  for 
teleseismically  determined  centroid  moment  tens(}rs  (CMT)  reported  by  Harvard  fcx  38  earthquakes 


39 


Figure  18.  Ccnnposite  focal  plot  of  teleseismically  determined  CMT  mechanisms  far 
earthquakes  which  may  occur  within  the  Efate  nest  Data  are  38  Harvard  CMT  occurring 
since  1977  with  epicenters  occurring  between  18.1*’$  aiul  17.0*’S  between  167 .0°E  and 
168.1*’E.  Circles  are  P  axes,  triangles  are  T  axes,  and  pluses  are  B  axes  for  individual 
CMT,  the  heavy  line  shows  the  composite  mechanism  formed  by  adding  together  the 
individual  moment  tensors  (see  text). 
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occurring  since  1977  in  this  aietu  These  CMT  are  remarkable  in  that  all  but  three  are  almost 
idbitical  to  one  another.  For  the  38  CMT  the  seismic  consistency  Cs  as  de^ed  by  Frohlich  and 
Apperson  (1992)  is  0.97,  i.e..  if  we  form  the  composite  tensor  by  adding  together  the  38 
individual  CMT,  the  scalar  moment  of  the  composite  tensor  divided  by  the  sum  of  individual  scalar 
moments  is  0.97.  Fm  the  ccHnposite  tensor  the  P  axis  has  a  dip  of  18**,  the  T  axis  has  a  dip  of 
about  72**  and  die  B  axis  is  horizontal  with  a  strike  of  14**W  of  N;  thus,  the  nodal  planes  have  dips 
ci  63**  westward  and  27**  eastward.  Neither  of  these  corresponds  to  the  possible  westward-dipping 
structure  with  a  dip  of  33**  observable  in  Figure  17. 
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